Intrinsic Coupling between Current and Domain Wall Motion in (Ga,Mn)As 
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We consider current-induced domain wall motion and, the reciprocal process, moving domain wall- 
induced current. The associated Onsager coefficients are expressed in terms of scattering matrices. 
Uncommonly, in (Ga,Mn)As, the effective Gilbert damping coefficient a w and the effective out-of- 
plane spin transfer torque parameter f3 m are dominated by spin-orbit interaction in combination with 
scattering off the domain wall, and not scattering off extrinsic impurities. Numerical calculations 
give a w ~ O.Of and f3 w ~ 1 in dirty (Ga,Mn)As. The extraordinary large j3 w parameter allows 
experimental detection of current or voltage induced by domain wall motion in (Ga,Mn)As. 



The principle of giant magneto resistance is used to 
detect magnetic information. Large currents in mag- 
netic nanostructures can manipulate the magnetization 
via spin transfer torques [l[ . A deeper knowledge of the 
coupled out-of-cquilibrium quasi-particle and magnetiza- 
tion dynamics is needed to precisely control and utilize 
current-induced spin transfer torques. 

The magnetization relaxes towards its equilibrium con- 
figuration by releasing magnetic moments and energy 
into reservoirs. This friction process is usually described 
by the Gilbert damping constant a in the Landau- 
Lifshitz- Gilbert (LLG) equation. Spins traversing a mag- 
netic domain wall exert an in-plane and an out-of-plane 
torque on the wall [2|]. In dirty systems, when the domain 
wall is wider than the mean-free-path, the out-of-plane 
torque, often denoted the non-adiabatic torque, is pa- 
rameterized by the so-called /3-factor 0]. The Gilbert 
damping coefficient a, the in-plane spin-transfer torque, 
and the out-of-plane torque coefficient (3 determine how 
the magnetization is influenced by an applied current, 
e.g. the current-induced Walker domain wall drift veloc- 
ity is proportional to (3 /a @, S 0] • Scattering off impu- 
rities are important for a and (3 Additionally, 
domain wall scattering can contribute to a and (3. In 
ballistic (Ga,Mn)As, intrinsic spin-orbit coupling causes 
significant hole reflection at the domain wall, even in the 
adiabatic limit when the wall is much thicker than the 
Fermi wavelength [5|. This grossly increases the out-of- 
plane spin-transfer torque, and consequently the current- 
driven domain wall mobility. So far, there are no inves- 
tigations on the effect of these domain wall induced hole 
reflections on the effective Gilbert damping constant a. 

Experimental (Ga,Mn)As samples are dirty so that the 
effect of disorder on the effective Gilbert damping and the 
out-of-plane spin transfer torque should be taken into 
account. We find surprisingly that, in systems with a 
large intrinsic spin-orbit coupling, domain wall scattering 
contributes dominantly to a and (3 even in the dirty limit. 
Intrinsic current-domain wall motion coupling is robust 
against impurity scattering. 

Current-induced domain-wall motion has been seen in 
many experiments Q. The reciprocal effect, domain- wall 
motion induced current, is currently theoretically investi- 



gated [1,0], and seen experimentally [||. A precessing do- 
main wall induces a charge current in ferromagnetic met- 
als 6] similar to spin-pumping in layered ferromagnet- 
normal metal systems [9( . For rigid domain wall motion, 
the induced charge current is proportional to (3 /a 0- We 
find that (3 and (3 /a in (Ga,Mn)As are so large that the 
current, or equivalently, the voltage induced by a moving 
domain wall is experimentally measurable. 

Onsager's reciprocity relations dictate that response 
coefficients of domain wall motion induced current and 
current induced domain wall motion are related. In dirty 
systems, these relations have been discussed in Ref. [7|. 
Ref. [7fl also used the scattering theory of adiabatic pump- 
ing to evaluate the non-adiabatic spin-transfer torque in 
ballistic systems without intrinsic spin-orbit interaction. 
We first extend the pumping approach to (Ga,Mn)As 
with strong intrinsic spin-orbit interaction, and second, 
also evaluate the Onsager coefficient as a function of 
sample disorder. In determining all Onsager coefficients, 
magnetization friction must be evaluated on the same 
footing. To this end, we generalize the energy pump- 
ing scattering theory of Gilbert damping [1 01 ] to domain 
wall motion. Our numerical calculation demonstrates, 
for the first time, that domain wall scattering is typi- 
cally more important than impurity scattering for the 
effective domain wall motion friction in systems with a 
strong intrinsic spin-orbit interaction. Our novel Onsager 
scattering approach can also be used to compute the ef- 
fective rigid domain wall motion a and (3 parameters in 
realistic materials like Fe, Ni, Co, and alloys thereof from 
first-principles. 

Let us discuss in more detail the Onsager reciprocity 
relations in our system. The magnetic field is a thermo- 
dynamic force for the magnetization since it can move 
domain walls. The electric field is a thermodynamic force 
for the charges as it induces currents. In systems where 
charge carriers also carry spin, the magnetic and charge 
systems are coupled. Through this coupling, the elec- 
tric field can move a domain wall and, vice versa, the 
magnetic field can induce a current. This phenomenon, 
where the thermodynamic force of one system can induce 
a flux in another system is well-known in thermodynam- 
ics ll|: Assume a system described by the quantities 
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{qi}, Xi denotes the thermodynamic force, and Ji the 
flux associated with the quantity qi. In linear response, 
Ji = LijXj, where Ly are the Onsager coefficients. 
Onsager's reciprocity principle dictates — CiCjLji, 
where u = 1 (e^ = —1) if is even (odd) under time- 
reversal (llj . Fluxes and forces are not uniquely defined, 
but the Onsager reciprocity relations are valid when the 
entropy generation is S = J2i JiXi [H|. 

We first derive expressions for the Onsager coefficients 
and determine the Onsager reciprocity relations between 
a charge current and a moving domain wall in terms 
of the scattering matrix. Subsequently, we derive the 
relation between the Onsager coefficients and the effec- 
tive Gilbert damping parameter a w and the out-of-plane 
torque parameter j3 w for domain wall motion. Finally, 
we numerically compute a w and (3 W for (Ga,Mn)As. 

We start the derivation of the Onsager coefficients in 
terms of the scattering matrix by assuming the following 
free energy functional for the magnetic system 

F[M] = M s Jdr(J^[(V6) 2 + sin 2 (0)(V<£) 2 ] + 

^ sin 2 (0) sin 2 (0) - ^ cos 2 (0) - H ext cos(0) J , (1) 

where M s , J and H QXt are the saturation magnetization, 
spin-stiffness and external magnetic field, respectively, 
and K z and K± are magnetic anisotropy constants. The 
local magnetization angles 9 and cf> are defined with re- 
spect to the z- and x-axis, respectively. The system con- 
tains a Bloch wall rotating in the (transverse) x-z plane, 
008(6*) = ta,uh([y-r w ]/X w ), sin(6>) = 1/ cosh.([y-r w ]/X w ), 
where r w is the position of the wall, and A^, is the wall 
width. We assume the external magnetic field is lower 
than the Walker threshold, so that the wall rigidly moves 
(<f> = 0) with a constant drift velocity. In this case r w 
and <j) completely characterize the magnetic system, and 

A,„ = 



The current is along 



\/j/(K z 
the y-axis. 

The heat dissipated per unit time from a charge current 
J c is Q = J c (Vl — Vr), where Vl (Vr) is the voltage in 
the left (right) reservoir. Using the relation dS = dQ/T, 
this implies an entropy generation S — J c (Vl — Vr)/T. 
Thus, X c = (Vl — Vr)/T is the thermodynamic force 
inducing the flux J c . We assume the magnetic system 
to be at constant temperature, which means that the 
heat transported out of the magnetic system as the do- 
main wall moves equals the loss of free energy. This 
implies an entropy generation S = Q/T = —F/T — 
(—dF[r w , 4>]/Tdr w ) r w — X W J W , where we have defined 
the force X w = ~dF[r w , 4>]/Tdr w and flux J w = r w . Us- 
ing Eq. Q, we find X w = -2AM s H cxt /T, where A is the 
conductor's cross-section. Fluxes are related to forces by 



Jc 



(2) 
(3) 



where L cc — GT and G is the conductance. L ww (L wc ) 
determine the induced domain wall velocity by an exter- 
nal magnetic field (a current). The induced current by a 
moving domain wall caused by an external magnetic field 
H cxt is controlled by L cw . Both charge and r w are even 
under time-reversal so that L cw — L wc 

The current induced by a moving domain wall is para- 
metric pumping in terms of the scattering matrix [9j : 
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where S a g is the scattering matrix between transverse 
modes in lead /3 to transverse modes in lead a. The 
system has two leads (a,/3 G {1,2}). The trace is over 
all propagating modes at the Fermi energy Fp. From 
Eqs. @ and ([3]) we find J c = L cw f w /L ww . 

We consider transport well below the critical transition 
temperature in (Ga,Mn)As, which is relatively low, and 
assume the energy loss in the magnetic system is trans- 
ferred into the leads by holes. Generalizing Ref. [T(| to 
domain wall motion, this energy-flux is related to the 
scattering matrix: 
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For a domain wall moved by an external magnetic field, 
we then find that X W J W = J%,/L ww — Je/T. In sum- 
mary, the Onsager coefficients in Eq. ^ and Eq. are 
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where t is the transmission coefficient in the scatter- 
ing matrix. We have omitted the temperature factor 
in the coefficients ([6]), ([7]), and ([8]) since it cancels with 
the temperature factor in the forces, i.e. we transform 
L — > L/T and X — * TX. The Onsager coefficient ex- 
pressions in terms of the scattering matrix are valid irre- 
spective of impurity disorder and spin-orbit interaction in 
the band structures or during scattering events, and can 
treat transport both in ballistic and diffusive regimes. 

Let us compare the global Onsager cofficients ((6|), (JTJ, 
and (jl]) with the local Onsager coefficients in the dirty 
limit to gain additional understanding. In the dirty limit, 
all Onsager cofficients become local and the magnetiza- 
tion dynamics can be described by the following phe- 
nomenological local LLG equation [2, 0] : 



-7m x H g + am x m 
- (1 — /3mx) (v s • V) m, 



(9) 
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where m is the magnetization direction, H e ff is the 
effective magnetic field, 7 is the gyromagnetic ratio, 
v s = —~hPj/(eS ), Sq — M a /7, M s the magnetization, 
a the Gilbert damping constant, P the spin-polarization 
along — m of the charge carriers and (3 is the out- 
of-plane spin-transfer torque parameter. Substituting 
a Walker ansatz into Eq. ([9]) gives below the Walker 
threshold ji]: ar w /X w = -7-ffext - iWj/ (eSbA w ). In 
dirty, local, systems this equation determines the rela- 
tion between the flux J w and the forces X w and X c 
as L ww = \ w /(2ASqo) and L wc — —h/3PG / (eaSoA) , 
where we have used j — <t(Vl — Vr)/L, and G = aA/L. 
Here, L is the length of the conductor, e the electron 
charge, and a the conductivity. This motivates defining 
the following dimensionless global coefficients: 

^ X w c L wc 

1ASqL ww 2hPG L ww 

a w is the effective Gilbert damping coefficient and (3 W is 
the effective out-of-plane torque on the domain wall. 

We will in the following investigate a w and (3 W for 
(Ga,Mn)As by calculating the scattering matrix expres- 
sions in Eq. (J6j) and Eq. ([7|) . We use the following Hamil- 
tonian to model quantum transport of itinerant holes: 

H = H L + h(r) • J + V(r) . (10) 

Here, £/l is the 4x4 Luttinger Hamiltonian (parame- 
terized by 71 and 72) for zincblende semiconductors in 
the spherical approximation, while h ■ J describes the 
exchange interaction between the itinerant holes and the 
local magnetic moment of the Mn dopants. We introduce 
Anderson impurities as V(r) = Y^i^i^Rt, where Ri is 
the position of impurity i, Vi its impurity strength, and 
5 the Kronecker delta. More details about the model and 
the numerical method used can be found in Refs. [Til, flit- 

We consider a discrete conductor with transverse di- 
mensions L x = 23nm, L z = Ylnm and length L y — 
400nm. The lattice constant is lnm, much less than the 
typical Fermi wavelength A^ ~ 8nm. The Fermi energy 
Ep = 82meV is measured from the bottom of the lowest 
subband. |h| = 0.5 x 10~ 20 J and 71 = 7. The typical 
mean-free path for the systems studied ranges from the 
diffusive to the ballistic regime / ~ 23nm — > 00, and we 
are in the metallic regime kpl ;§> 1. The domain wall 
length is X w — AOnm. The spin-density So from the local 
magnetic moments is So — l0hx/a,Q aAs , aQaAs the lattice 
constant for GaAs, and x = 0.05 the doping level [l4j|. 

Fig. shows the computed effective Gilbert damping 
coefficient a w versus X w /l for (Ga,Mn)As containing one 
Bloch wall. Note the relatively high a w ~ 5 x 10~ 3 in 
the ballistic limit. Additional impurities, in combination 
with the spin-orbit coupling, assist in releasing energy 
and angular momentum into the reservoirs and increase 
a w . However, as shown in Fig.[T^,, impurities contribute 
only about 20% to a w even when the domain wall is two 
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FIG. 1: (a): Effective Gilbert damping a w as function of 
X w /l, where X w is the domain wall length and I is the mean 
free path when 72 = 2.5. Here, A„, is kept fixed, and I is 
varied. Inset: a w as a function of spin-orbit coupling 72 for a 
clean system, I — 00. (b): /3 W as a function of X m /l, where X m 
is the domain wall length and I is the mean free path when 
72 = 2.5. Here, X w is kept fixed, and I is varied. Inset: f3 w as 
a function of spin-orbit coupling 72 for a clean system, I — 00. 
In all plots, line is guide to the eye. 



times longer than the mean free path. Due to the strong 
spin-orbit coupling, ballistic domain walls have a large 
intrinsic resistance [B| that survives the adiabatic limit. 
When itinerant holes scatter off the domain wall their 
momentum changes and through the spin-orbit coupling 
their spin also changes. This is the dominate process for 
releasing energy and magnetization into the reservoirs. 
The saturated value « ffi ~6x 10 -3 is of the same order 
as the estimates in Ref. [l6| for bulk (Ga,Mn)As. The 
inset in Fig. [T^, shows the domain- wall contribution to 
a w versus the spin-orbit coupling for a clean system with 
no impurities. a w monotonically decreases for decreasing 
72 and vanish for 72 — > 0. Since, X w /Xf ~ 5, itinerant 
holes will, without spin-orbit coupling, traverse the do- 



4 



main wall adiabatically. 

Fig. [TJd shows j3 w versus X w /l. (3 W decreases with 
increasing disorder strength. This somewhat counter 
intuitive result stem from the fact that domain walls 
in systems with spin-orbit coupling have a large intrin- 
sic domain wall resistance [B| which originates from the 
anisotropy in the distribution of conducting channels [B| . 
The reflected spins do not follow the magnetization of 
the domain wall, and thereby cause a large out-of-plane 
torque Q. This causes the large (3 W in the ballistic limit. 
Scalar, rotational symmetric impurities tend to reduce 
the anisotropy in the conducting channels, and thereby 
reduce the intrinsic domain wall resistance and conse- 
quently reduce (3 W . Deeper into the diffusive regime, j3 
saturates. Here, the domain wall resistance and (3 W are 
kept at high levels due to the increase in the spin-flip 
rate caused by impurity scattering. The saturated value 
is /3 ~ 1. For even dirtier systems than a reasonable 
computing time allows, we expect a further increase in 
(3 W . In comparison, simple microscopic theories for fer- 
romagnetic metals where one disregards the spin-orbit 
coupling in the band structure predict f3 ~ 0.001 — 0.01 
[1, [!, 3 . Similar to the Gilbert damping, in ballistic sys- 
tems f3 w increases with spin-orbit coupling because of the 
increased domain wall scattering i, see Fig. [lb inset . 

/3 W can be measured experimentally by the induced 
current or voltage from a domain wall moved by an ex- 
ternal magnetic field as a function of the domain wall 
velocity [7[. From the Onsager relations we have that 
J c = L CW X W . Using X w — J w /L ww , the induced current 
and voltage are Q: 

h PC hP 
J c = -2/3— — r w => V = -20 w —f w . (11) 

An estimate of the maximum velocity of a domain wall 
moved by an external magnetic field below the Walker 
treshold is r w ~ 10 m/s [IT]]. With X w = 40 nm and 
P = 0.66 this indicates an experimentally measurable 
voltage V ~ 0.2 fxV. 

In conclusion, we have derived Onsager coefficients and 
reciprocity relations between current and domain wall 
motion in terms of scattering matrices. In (Ga,Mn)As, 
we find the effective Gilbert damping constant a w ^0.0l 
and out-of-plane spin transfer torque parameter [3 W ~ 1 . 
In contrast to ferromagnetic metals, the main contribu- 
tions to a w and (3 W in (Ga,Mn)As are intrinsic, and in- 
duced by scattering off the domain wall, while impurity 
scattering is less important. The large (3 W parameter im- 
plies a measurable moving domain wall induced voltage. 
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